A topologically flat thick 2-brane on higher dimensional black hole backgrounds 
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We present a numerical solution for a topologically flat 2-dimensional thick brane on a higher 
dimensional, spherically symmetric black hole background. Present solution is the last, missing part 
of the complete set of solutions for the thickness corrected brane - black hole problem in arbitrary 
number of dimensions. We show that the 2-dimensional case is special compared to all the higher 
dimensional solutions in the topologically Minkowskian family as being non-analytic at the axis of 
the system. We provide the numerical solution in the near horizon region and make a comparison 
with the infinitely thin case. 
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I. INTRODUCTION 

The study of higher dimensional black holes, branes 
and their interactions is an active field of research in sev- 
eral different areas of modern theoretical physics (THl]- 
One interesting direction, which has been first introduced 
by Frolov 0, is to consider a brane - black hole (BBH) 
toy model for studying merger and topology changing 
transitions in higher dimensional classical general rela- 
tivity 0, [f| , or in certain strongly coupled gauge theories 
[3, 0] through the AdS/CFT correspondence [8|. Gener- 
alizations of the BBH model by studying thickness cor- 
rections to the Dirac-Nambu-Goto effective brane action 
P~[Tl| from higher order curvature terms, have also been 
studied recently, first by perturbative approaches (l2l . [l3j . 
and later within an exact description [14j |. 

The results of the perturbative approaches concluded 
that there is a "symmetry breaking" between the two 
topologically different solution family, as regular pertur- 
bative solutions do not exist for Minkowskian embed- 
ding topologies except in the special case of a 2-branc. 
The problem, on the other hand, can be solved regu- 
larly for any brane dimensions in the black hole embed- 
dings. This virtual "symmetry breaking" phenomenon 
obtained a simple resolution in fl3j , where it was pointed 
out that perturbative thick solutions break down around 
their thin counterparts for Minkowski topologies, because 
the thin solutions are not analytic at the axis of the 
system. Motivated by this observation, in 14 \, a gen- 
eral family of thick solutions could be provided for both 
topologies within a non-pcrturbativc numerical approach 
for all but one exceptional case. The exception, mysteri- 
ously, turned out to be the same case, where the regular 
perturbative solution existed, namely the 2-dimensional, 
topologically flat case. 

The above findings of [l4[ naturally raised the ques- 
tion: How can a regular perturbative solution exist in 



the same single case where the regular non-perturbative 
solution can not be found? 

In the present paper we provide the answer to this 
question and obtain the so far missing solution of the 
topologically flat 2-brane in arbitrary number of bulk di- 
mensions for the thick-BBH system. By including this 
solution, we complete the full set of solutions of the prob- 
lem that we started presenting first with a perturbative 
approach in |13| and continued with a non-perturbative 
description in [l4| . The present work, therefore, is the 
final part of the series of papers we addressed to the 
thickness corrected BBH problem and we kindly refer 
the reader also to [l3|, [l4[ for the more detailed model 
setups and for all those definitions, notations and results 
that might be missing here and would make the present 
paper completely self-contained. 

The plan of the paper is as follows. In Sec. II, we pro- 
vide a short overview on the thin- and thick-BBH model 
setups. In Sec. Ill we obtain the special case of the 2- 
brane equation, and in Sec. IV we analyze its regularity 
conditions for flat topology. In Sec. V we provide the 
non-perturbative, numerical solution of the problem in 
the near horizon region, and in Sec. VI we draw our con- 
clusions. 



II. THE BRANE - BLACK HOLE SYSTEM 

Let us overview quickly the most important properties 
of the thin-BBH system introduced by Frolov in H, and 
its thickness corrected generalization provided in [13[ . 

A. The thin model 

We consider static brane configurations in the back- 
ground of a static, spherically symmetric bulk black hole. 
The metric of an ./V-dimensional, spherically symmetric 
black hole spacetime is 
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ds 2 



9ab 



dx a dx b 



fdt 2 + f- dr 



(l) 
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where / = f(r) and <Kl 2 N _ 2 is the metric of an N — 
2 dimensional unit sphere. One can define coordinates 
9i(i = 1, . . . , N — 2) on this sphere with the relation 



dn 2 l+1 = d6 2 i+1 + sin 2 9 i+1 dQ 2 



(2) 



The explicit form of / is not important, it is only assumed 
that / is zero at the horizon ro, and it grows monotoni- 
cally to 1 at the spatial infinity r —¥ oo, where it has the 
asymptotic form (l5j . 



/ = 1 



r.JV-3 



(3) 



In the zero thickness case the test brane configura- 
tions in an external gravitational field can be obtained by 
solving the equation of motion coming from the Dirac- 
Nambu-Goto action (9T-lllj|. 



S = 



d D ' 



C A /-det7 



\JLV ' 



where j„ v is the induced metric on the brane 
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fJ,V 



9ab 



dx a dx b 
WW 



(4) 



(5) 



and C m (a* — Oj ■ • ■ 7 D — 1) are coordinates on the brane 
world sheet. The brane tension does not enter into the 
brane equations, thus for simplicity it can be put equal 
to 1. It is also assumed that the brane is static and 
spherically symmetric, and its surface is chosen to obey 
the equations 



On-2 



ir/2 



(6) 



With the above symmetry properties the brane world 
sheet can be defined by the function 6d—i = 9{r) and we 
shall use coordinates C M on the brane as 



C = {t, 



. . , cj) n } with n = D — 2 



(7) 



The parameter n denotes the number of dimensions in 
which the brane is rotationally symmetric. In this paper 
we consider the special case of n = 1, i.e. a 3-dimcnsional 
brane world sheet, that is a 2-dimensional, axisymmetric 
brane embedded into the higher dimensional black hole 
spacetime. 

With this parametrization the induced metric on the 
brane is 



j^d^dC = -fdt 2 



dr 2 +r 2 sin 2 (Ml 2 , 



(8) 

where, and throughout this paper, a dot denotes the 
derivative with respect to r, and the action (j4|) reduces 
to 



S = MA n Co dr , 



(9) 
(10) 

where At is the interval of time and A n = 2n n / 2 /T(n/2) 
is the surface area of a unit n-dimensional sphere. 



C = r n sin™ 0yf 1 + fr 



B. Thickness corrections 

In the case of a thick brane, the curvature corrected 
effective brane action is obtained by Carter and Gregory 
in [ill, and the corrections to the thin DNG action are 
induced by small thickness perturbations as 



S = I d D (J-det>y 



C X R + C 2 K 2 



(11) 



where R is the Ricci scalar, K is the extrinsic curvature 
scalar of the brane and the coefficients C\ and Ci are 
expressed by the wall thickness parameter t as 



C^^i 2 



Co 



1 



24 3 
The parameter fi is related to the thickness by 

1 



(12) 



(13) 



which originates from a field theoretical domain-wall 
model where [i is the mass parameter and A is the cou- 
pling constant of the scalar field. 

After integrating out the spherical symmetric part and 
the time dependence on the introduced static, spherically 
symmetric, higher dimensional black hole background, 
one obtains (see also [l3T |) 



(14) 
(15) 



S = MAn I 

8/i 2 



C 



'if 



■C [l + eS] , 



where we introduced the notations 

i>2 



L 2 



S = aK z + bQ , 



with 



Q = K a h K h a , a 



n 2 -U 
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-L , b 



6-TT 2 



24 



-L 



(16) 



(17) 



Here L is the relevant dynamical length scale of the sys- 
tem which has to be large compared to the thickness 
parameter £ in order to (|11|) remain valid. The explicit 
expressions of the curvature scalars K and Q are given 
in (35) and (36) of 0. 

For a detailed introduction of both the thin- and thick- 
BBH systems, please refer to (2T [l3j. 



III. THE 2-BRANE EQUATION 

From this section on we will focus only on the case of 
the topologically flat or Minkowskian 2-brane. In two di- 
mensions "topologically flat" is synonymous with "topo- 
logically Minkowskian" , and we retain both terms for the 
sake of elegant variation. 
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In order to obtain the 2-brane equation of motion, first 
we observe that the thickness corrected DNG-brane ac- 
tion is a function of the second derivative of 9 and thus 
the Euler-Lagrange equation of the problem has the form 
(see for example jl7j |) 



d 2 




- d 1 


fdC\ 


dC 


dr 2 


\d9j 


dr 1 


Kd9j 


' 39 



= 



From (fT5|) the actual equation of motion becomes 

0(4) + Tl0 (3) + ^ ^ ^ ^ /(3 ) ^ £ /, /, r ) = , 

where, in the 2-dimensional (n = 1) case 



(18) 



(19) 



with 



[ 6/ + 4r/ + 2rf cot I 
2r 3 / 2 cot 99 3 - 10r 3 f 2 90} 



rfF 2 



r 2 f Uf + rf 



F=Jl + fr 2 6, 



(20) 



(21) 



and Ti is given in the Appendix. 

As it is immediate to sec, (|19p is a 4i/i-order, highly 
nonlinear equation, and it is probably impossible to 
present its solutions in closed, analytic form. Hence, the 
goal of this paper, is to provide a regular, numerical solu- 
tion of (flU)) in arbitrary bulk dimensions for flat topology. 



IV. REGULARITY AND BOUNDARY 
CONDITIONS 

It was pointed out first in Il3j, that the brane solu- 
tions obtained by Frolov in [2| are regular but not an- 
alytic (or smooth) at the axis of the thin-BBH system 
for the Minkowski embedding branch. In fact they are 
not even differentiable at that point and thus belong to 
the class of C° functions only. According to this prop- 
erty, and since we found that the perturbative approach 
broke down around these solutions for the thick case, we 
concluded that the thick brane solutions must behave 
significantly differently at the axis of the system, namely 
we expected them to be smooth there. It was surprising 
however that in the single case of the 2-brane, a regular 
perturbative solution existed. We gave a detailed analy- 
sis of this solution in [l]| . 

In [Tij . approaching the problem by a new, non- 
perturbative numerical method, we looked for the miss- 
ing solutions of the Minkowski branch in the class of ana- 
lytic functions. We obtained regular boundary conditions 
at the axis of the system by considering the series expan- 
sion of the exact 4i/i-order equation of motion around 
9 = 0. With this method we successfully provided all 
the missing, topologically Minkowskian solutions of the 
thick-BBH system, except in the curious case of the 2- 
brane again. 

It was obvious, of course, that a perturbative solution 
can not exist without the existence of a non-perturbative 



solution, and also since explicitly constructed, field theo- 
retical domain wall solutions [l8|, [l9| clearly exist in the 
case of the 2-brane, we suspected that the lack of this so- 
lution must lie somewhere in the validity of the applied 
method. Nevertheless, we were so enthusiastic with pro- 
viding the whole family of the missing solutions in the 
C°° class at the axis, that it didn't occur in our mind at 
the time, that the 2-dimensional case might be special in 
the topologically Minkowskian family as being the only 
one which is non-analytic at the axis. 

The main result of the present paper is the observation, 
that the 2-dimensional solution of the thick-BBH system 
is in fact a special one in the topologically Minkowskian 
family as being C° (or as we will soon see maximum C 1 ) 
function at the axis of the system. All the other dimen- 
sional solutions arc C°° . This also explains the fact why 
it was only the 2-dimensional case where a perturbative 
solution could exist around the thin solution. 

In the remaining of this section we analyze the asymp- 
totic behavior of (|T9"t near the axis of the system. We 
obtain necessary boundary conditions from regularity re- 
quirements and show that these conditions can always be 
fulfilled in order to obtain a regular solution. 



Asymptotic analysis 

We know from the above considerations that a regu- 
lar solution of the problem must exist although analytic 
solution could not be found at the axis of the system. 
Thus the point r\ on the axis, where 9(r{) = 0, must be 
a regular singular point of the differential equation (|19p . 
Even though Q19[) is highly nonlinear, general results from 
the theory of local analysis of linear differential equations 
can be applied, because we know from physical considera- 
tions that (|19[) should not develop any nontrivial singular 
points in its domain. 

Hence, if a solution is not analytic at a regular singular 
point (see e.g. [13), its singularity must be either a pole 
or an algebraic or logarithmic branch point, and there is 
always at least one solution of the form 



(r) = (r - ri ) a A(r) 



(22) 



where a is a number called indical exponent and A(r) is 
a function which is analytic at n and has a convergent 
Taylor series. 

In the general case a can be any number that solves 
(|19[). In our specific case however, the thin solution, 
around which a regular perturbative solution existed, had 
the asymptotic form near n (see 0, EH ) 



9(r) = r)\Jr - n + 



(23) 



and similarly, the corresponding perturbative solution 
(see (HI) near the same point had the asymptotic form 



hhick = 9 t hin +eip= (rj + £K)\/r-ri + ..., (24) 
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where (p is the perturbation function and r\ and k are 
coefficient functions defined in (T^. Thus, in order to 
obtain the non-perturbative thick solution, we also chose 
a to be 1/2, in accordance with the perturbative results. 
This choice will also have the advantage of naturally fix- 
ing the free boundary condition in the next section for a 
unique numerical solution. 

With a = 1/2, from (|2"2"j) we obtain that the asymptotic 
form of 9(r) near the axis is 



(r) = Axy/r - n + A 2 (r - n) ' 
+ A 3 (r - n)5 + A 4 (r - n)i + 



(25) 



Plugging this expression and its derivatives into (|T9|) . one 
can obtain the following asymptotic behavior 



C3 



Cl 



(r-n)5 (r-n)= \[r - ri 



c 



= 0, (26) 



where the coefficient functions Ci are polynomial expres- 
sions of Ai with dependences 

c 5 = c 5 (A 1 ,A 2 ) 

c 3 = c 3 (Ai,A 2 ,A 3 ) 

ci = ci(Ai,A 2 , A 3 , Aa) . 

In order to obtain a regular solution, we need to require 
for the Ci coefficient functions to disappear at n. From 
the explicit forms of Ci, one can find that the coefficient 
A\ can be chosen freely, and once it's fixed, the remain- 
ing coefficients A 2 , A3 and A4 can be computed from the 
requirements that c^(ri) = 0. Consequently, the solu- 
tion is uniquely determined by the parameter r\ (i.e. the 
minimal distance parameter), and the explicit values of 
the coefficients A2, A3 and A4 can be immediately ob- 
tained from the c^(ri) = equations. (This procedure is 
of course necessary before the numerical setup.) For the 
question of existence, the Ci{r\) = equations arc always 
soluble, as it turns out that these are linear equations for 
the Ai coefficients. The explicit forms of the coefficients 
A2, A3 and A4, as successive functions of A\ and r\ are 
given in the Appendix. 

In conclusion we found that a regular 2-dimensional so- 
lution can always be given for the exact problem, and a 
unique solution is completely determined by the regular- 
ity requirements once the coefficient A\ is fixed. Since we 
are free to choose A\, it can also put to be for example. 
In this case the asymptotic form (|23|) starts with 



") = A 2 (r - n)5 + 



(27) 



which solution is a C 1 function, however in all other cases 
the solution of (HU) is C° at r x . 



;ained from the analysis of the thick- 
141, it is not too difficult to obtain 



the experiences we 
BBH system in 
the numerical solution here, after the initial conditions 
for this specific case has been clarified. 

As we discussed earlier, we are free to choose the coef- 
ficient Ai in the asymptotic solution (f2"5")l . Nevertheless, 
in order to be completely consistent with our previous 
perturbative results in the flat 2-brane case, we make the 
choice 



A\ = T] + EK 



(28) 



that was forced upon us by regularity requirements for 
the perturbations. Having fixed this freedom, the re- 
maining three conditions A 2 , A3 and A4 are determined, 
as discussed in the previous section, and the correspond- 
ing numerical solution is unique. 

In obtaining the numerical solution we used 
Mathcmatica® NDSolve function. The integration 
range went from r% until 1000 to check the accordance 
with the corresponding perturbative solution. 

The configurations of the perturbative thick 2-brane 
solutions in the near horizon region have been presented 
in [l]|. Since the effects of the nonlinearities are really 
small even in the gravitationally strong, near horizon re- 
gion, the overall global picture of the thickness corrected 
brane configurations remain very similar to the pertur- 
bative case. To make however the small effects visible, 
we plot on Fig. 1 and Fig. 2 the difference function 

A6(r) = 6(r) - 6 DNG (r), 

of the present thick- and the original thin solutions in 
4 and 5 dimensions. A8(r) is the exact analog of the 
perturbation function eip(r) defined in (f24)). 



A6 




FIG. 1. The picture shows a sequence of near horizon A9(r) 
curves in 4 dimensions with minimum horizon distance range 
1.01 < n < 2. 



V. NUMERICAL SOLUTION NEAR THE 
HORIZON 

For illustrating the obtained results, we provide the nu- 
merical solution of the 2-dimensional flat problem. With 



Comparing the results with the corresponding plots 
of the perturbative solutions in [l3| . we find that the 
general behavior of the A8(r) curves are essentially the 
same. This is of course what one expects. On the other 
hand one also expects sonic differences coming from the 
nonlinear regime of (|19[) . and those are also present if 
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AO 




FIG. 2. The same sequence of A8(r) curves as on FIG.l in 5 
dimensions. 



wc enlarge the very near horizon region of the individual 
curves. The main features of these differences arc plotted 
on Fig. 3 and Fig. 4. 



in the perturbative case, the A6(r) curve has an extra lo- 
cal pattern near the minimum horizon distance r\. The 
extra pattern is that the curve goes through some addi- 
tional local extrema in this region before it finally tends 
to reach the global minimum which is similar to the one 
in the perturbative case. During these local differences 
however, A8(r) remains negative in this whole region. 

On the qualitatively similar curve on Fig. 4, the essen- 
tial difference compared to Fig. 3 is that A9(r) develops 
a sign change in this region, that is there are several in- 
tersections in the very near horizon region between the 
thick and thin solutions compared to the single inter- 
section that is present in the case of the perturbative 
solution (see Fig. 6 in fl3|). 

With increasing horizon distance r, the above nonlin- 
earities decay quickly, and the solution agrees very well 
with the perturbative results. 



VI. CONCLUSIONS 
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FIG. 3. Near horizon nonlinear effects on a A9(r) curve in 4 
dimensions with minimum horizon distance r\ = 1.19744. 
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FIG. 4. Near horizon nonlinear effects on a A9(r) curve in 4 
dimensions with minimum horizon distance ri = 1.2538. 

On Fig. 3 an interesting nonlinear effect appears com- 
pared to the perturbative solution. Instead of monoton- 
ically decreasing down to its global minimum, as (p does 



In the present work we studied the problem of a topo- 
logically flat 2-brane in a higher dimensional, thickness 
corrected BBH system. We provided a regular, non- 
perturbative, numerical solution for this special case 
based on earlier perturbative considerations [13j . The 
main result of this paper is the observation that the 2- 
dimensional case of the problem is a special one as being 
non-analytic at the axis of the system. This property 
makes it unique in the family of thick solutions, as in 
all other dimensions both the Minkowski and black hole 
embedding solutions are analytic in their entire domain. 

We analyzed the asymptotic behavior of the solution 
and obtained that it is at most a C 1 , but in the general 
case it is only a C° function at the axis, just like the 
corresponding thin solutions. The initial conditions of 
the problem are not uniquely fixed by regularity require- 
ments, thus the solution we provided here is not unique. 
It is however perfectly consistent with our earlier pertur- 
bative results. 

With the present paper, we have provided the complete 
set of solutions of the thick-BBH problem in a series of 
three consecutive papers. First, in [ill ], we obtained all 
possible perturbative solutions and later, in [l4| all non- 
perturbative solutions were given except the case of the 
flat 2-branc. Present work completes the set. 

In [l]| we analyzed the properties of a topology chang- 
ing, quasi-static phase transition in the thick-BBH sys- 
tem. The obtained results in the present case however 
does not change our previous findings and thus we don't 
consider the phase transition in this paper. 

The result, that thickness corrections change the an- 
alytic properties of the brane solutions at the axis of 
the system might have some physical consequences. In- 
finitely thin brane solutions, naturally, are very impor- 
tant in higher dimensional physics, but considering the 
present problem, one has the intuition that the thickness 
corrections, which are in agreement with field theoretic 
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domain wall models, made the thin-BBH system more 
stable in the sense that regular, analytic solutions could 
be provided in essentially all cases. Since small physi- 
cal perturbations to any system are usually proportional 
to the derivatives of the unperturbed solution, it is very 
possible that the thin-BBH solutions are not entirely sta- 
ble against small perturbations in the Minkowski branch. 
(Stability properties of the BBH system for the analytic 
black hole embedding solutions has been studied for ex- 
ample in (2lj.) This property however has been cured by 
the thickness corrections and it is somehow in accordance 
with our physical expectations. 

The special case of the 2-brane as "remaining" non- 
analytic after thickness corrections, thus, is an unex- 
pected property, which makes it physically interesting. 
So much the more that thick 2-branes, i.e. thin walls on 



black hole backgrounds in standard 4-dimcnsional gen- 
eral relativity are certainly real, physical objects. The 
fact that these solutions arc essentially different from the 
corresponding ones in higher dimensions is remarkable. 
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Appendix: Coefficient Functions 
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- (2ae - 26e + r 2 ) cos(26») + 8ae/r 2 sin 2 9)) - l/2/ 3 r 3 (-600(a + b)er 2 99 2 - 40(a + b)er 3 9 3 

+ 2e9 3 (-3{37a + 296) + 2(25a + 16b)r 2 9 cot 9) - 2r9*((-13a - b)e cot 9 + r 2 9(16ae + 16be + 3r 2 

+ 6(3a + 26)e/r 2 + 3(a + 6)ecot 2 9)) + 2r 2 5 (-35a£ - 236e - 15r 2 + (11a + 56)e/r 2 + 2(4a + 6)e/ (3 V 3 

- 6ae cot 2 9 - 126e cot 2 9 + 3ae esc 2 9 + 9be esc 2 9) - 2er9 2 '9 esc 6l(36(a + 6)r 2 6> cos 9 

+ (295a + 2716) sin 9) - r 3 9 & cot 9 esc 2 9(9ae + 21be - 4r 2 + (5ae - 76e + 4r 2 ) cos(26») 

- 2(17a + 46)e/r 2 sin 2 61)))) ] , 
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Aa = 188697600 ^ 5(a + fe)£r6 [ Z2A\r\- 15263640^ (a + b)er 3 + 1260A 2 er(-491aA 2 - 608A 2 6 (A.2) 

+ 2573aA 3 r + 2540A 3 6r) - 504(L4i A 2 er 2 (A 2 (497a + 5036) - 4606A 3 (a + b)r) 
+ 420Alr(-65ae - 386e + 6r 2 ) + Atr 3 {~6ae - 346e + 7r 2 ) + 14^r 2 (46ae + 406e + 15r 2 ) 
+ 21A b 1 A 2 r i (222ae - 38be + 65r 2 ) - 630yl?(2(71a + 176b)e + A 2 2 r 3 (27ae - 416e + 17r 2 )) 
+ 210 J 4 4 r 2 (v4 3 r(983ae + 10676e - 21r 2 ) + A 2 (490a£ + 1606e + 93r 2 ))) 

+ (l/(/ 3 ))7(-33177600(a + b)e + Air(9(Mi£/r(-512(616a + 6016) - 128A 2 (611a + 5816)/r 2 

- 32A 4 (67a + 1196)/V + 40,4?(3a + 6)/V - A? (a + 6)/V) + 15/(-1536(4Ai(646a + 6316)e 
+ 7383A 2 (a + b)er + A\r{l$(a + b)e + r 2 )) + A 2 r 2 (768Ai(267a + 2486)e/r 

+ 6Al(a + 6)e/ 4 r 6 (4yl 1 + 15A 2 r) + A 4 e/V(-48Ai(47a + 176) + 8(5a(A? - 96A 2 ) 
+ {A\ - 150A 2 )6)r + 2L4?(a + b)fr 3 ) + 16/(120A?r 3 + 3A?(491a + 6856)e/r 3 

- 8e(9138aAi + 8802Ai6 + 253aA?r + 28668aA 2 r + 163A?6r + 28038A 2 6r)) 
+ 8 J 4 2 / 2 r 2 (240 J 4i(9a - 56)e - 48A 2 (137a + 1486)er - 3,4? (83a + 29b)efr 3 

- 4A?r(53ae + 55be + 9r 2 )))) + 4A'i/ : V 4 (90A 2 e/ 2 r 2 (A 2 (41a + 156) + 15yliA 2 (25a + 96)r 
+ 5(18A 2 + 5A 1 A 3 )(3a + 6)r 2 ) + 3/(16A?(2o + 6)er 2 - 24480A 2 (2a + 6)er 2 

+ 720Aier(87aA 2 + 36A 2 6 - 14A 3 (2a + b)r) + 60A 4 r(5ae + 26e + 9r 2 ) + 30A 3 ir 2 (25A 3 r(2ae + r 2 ) 
+ 4A 2 (97ae + 416e + 30r 2 )) + 180,4 2 (-8(2a + b)e + 15Al(2aer 3 + r 5 )) + 15A?er 3 (2/(3A?(20a + 76) 
+ 30A 1 A 2 {7a + 36)r + 5(18A 2 + 5AiA 3 )(a + 6)r 2 ) + Air(Ai(a + 6)/ (4) r + 2/ (3) (Ai(lla + 56) 
+ 15A 2 (a + 6)r)))) + r(90A?e/ 2 r 3 (A 1 (9a + 56) + 15A 2 (a + b)r) + /(135A 4 (a + 6)e/ (3) r 4 
+ 270A 3 > 3 (2Ai + 5A 2 r) + 4e(1080A 2 (3a + 6) + 15^i(52aA? - 1788aA 2 + 15/^6 - 828A 2 6)r 
+ a(8Al + 2235,45% - 12240v4j - 5040^1 A 3 )r 2 )) + 15Air(-192aAie/ (4) r 

+ / (3) (-48Ai (19a + 6b)e - 2208ayl 2 er + yl?(38aer + 3r 3 ))))) - Ai/ 2 r(384(270A 2 (249a + 2396)e 
+ 619200A 2 (a + 6)er 2 + 30Aier(15A 2 (697a + 6836) - 6682A 3 (a + b)r) + 30A 4 r(77ae + 556e - 3r 2 ) 
+ 15AlA 2 r 2 (247ae + 2036e + 7r 2 ) - A\r 2 (23ae + 36e + 10r 2 )) + A 2 r 2 (90A 4 (a + 6)e/ 3 r 4 (6A 2 + 90A 2 r 2 
+ 5Air(12yl 2 + 5A 3 r)) + A?e/V(8(-360A?(8a + 36) + 30Ai(5aA? - 939aA 2 + A 3 ± b - 339A 2 6)r 
+ (At + 195A 3 A 2 - 1530A 2 - 630AiA 3 )(5a + 6)r 2 ) + 135A?(a + 6)r 3 (yli/ (3 V + /(8Ai + 30A 2 r))) 
+ 240Air(15Ai(49a + 396)e/ (3 V + /(24Ai(52a - 156)e - 6A 2 (457a + 5876)er - 3A?(19a + 76)e/V 3 

- 2A\r(7&ae + 50be + 9r 2 ))) + 6/(-720A?r 3 (2yli + 9yl 2 r) - 16e(60A?(105a - 46) + 20A 1 {25aA 3 1 

- 330aA 2 + 6Atb - 261yl 2 6)r + {a(7Af + 1080yl?A 2 - 30060^1 + 38265^iyl 3 ) + 3(3,4? + 190yl?yl 2 

- 9780A 2 + 13085yliA 3 )6)r 2 ) + 5A 3 er 3 (-24Ai(19a + 76)/ (3) r + /(-48Ai(51a + 196) 

+ 8(5aA? - 588aA 2 + A\b - 204A 2 6)r + QA\(a + 6)/r 3 ))))) - 24Ai/V(4(2a + 6)e(-180A? 
+ 1620AiA 2 r + (2A\ + 165AlA 2 - 3060A 2 - 1260Ai^ 3 )r 2 - 75A 2 (18A 2 + 5Ai^ 3 )r 3 ) 

- 60A 2 r 3 (3A 2 + 90A 2 r 2 + 5Air(9A 2 + 5A 3 r)) - 15A 2 er 2 (2(8a + 36)/( J 4 2 + 90A 2 r 2 
+ 5A ir (6A 2 + 5A 3 r)) + r(2/(2A 2 (lla + 36) + 15AiA 2 (14a + 36)r 

+ 10a(18^ + 5Ai A 3 )r 2 ) + A 1 r(2aA 1 / (4) r + / (3) (3^6 + 20a(^! + 3A 2 r)))))))) } | rj , 

M = n. ggn4 / , m t 2 I (-14592(a + b)e + A 2 r(3e/r(-16(29a + 356) + 8A?(5a + 6)/r 2 (A.3) 
26880Ai (a + b)sjr A 

- A\(a + b)f 2 r 4 )24Aif 2 r(2Air 3 + 2{2a + 6)e(12 + A\r) + A 2 er 2 ((8a + 36)/ + 2a/r)) 
+ 4/(-8(6(31a + 286)e + A\r(lVae + 7be + 3r 2 )) + 3A?r 2 (e/r(-32a + A\(a + b)fr 2 ) 
+ /(-48(2a + b)e + A?(3a + b)efr 2 + A\r(2ae + r 2 )))))) ] | , 
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M = n4048M 3 (a + b)£r 4 I (15s(56576(a + b) + 384A 2 (74a + 696)/r 2 + 32^(37a + 576)/V (A.4) 

- 8At(10a + 3b)f 3 r 6 + A\{a + o)/V))// 2 + 12Q A\ f r\&A\r 3 {A 1 + 5A 2 r) 

- 2(2a + 6)e(12Ai + A\t - MA 2 r - 15A 2 ,4 2 r 2 ) + A 2 er 2 ((8a + 3b)f(2A 1 + 15A 2 r) 
+ r{2aA 1 f (3) r + f(UaAi + 3A±b + 30aA 2 r)))) + (l//)5 J 4ir(64(30Ai(101a + 95b)e 
+ 9372A 2 (a + b)er + A\r(7Qae + 67be - 9r 2 )) + A\r 2 {-\^A X (35a + 33o)e/r 

- 3Af(a + 6)e/V(4Ai + 15A 2 r) + A 2 1 ef 2 r 2 (24A 1 (37a + 136) - A(A\ - 42yl 2 )(5a + b)r 

- 9Al(a + b)fr 3 ) + 16/(12Ai(-30a + 136)e + 3A 2 (517a + 523b)er + 12A?(3a + b)sfr 3 
+ Alr{40ae + 34be + 9r 2 )))) + 4 J 4 2 r 2 (4(18(L4 2 (25a + 32&) £ + 18(L4i A 2 (173a + 172b)er 
+ 123480^2 (a + b)er 2 + mA\A 2 r 2 (37ae + 65fe - 7r 2 ) + 3(L4ir(30ae + 16te + 3r 2 ) 

+ Atr 2 (6ae - Ube + 5r 2 )) + 5A?r 2 (3A 2 e/V(Ai(25a + 9fe) + 15^ 2 (3a + b)r) 

+ 3/(8Ai(29a + 126)e - 336A 2 (2a + b)er + VbA%A 2 r 2 {2ae + r 2 ) + 4Alr(7ae + 3be + 2r 2 ) 

+ A 2 1 er 3 (A 1 (a + fe)/ (3) r + /(2Ai(7a + 3b) + 15A 2 (a + b)r))) + r(-144aAie/ (3) r 

+ /(-72Ai(7a + 3b)e - 672aA 2 er + 3A({a + b)efr 3 + A\(22aer + 3r 3 ))))) ] | ri . 
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